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introduced in one-loop amplitudes. Finally, we show how to derive transition 
rules between different schemes without having to do any loop calculations. 
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1 Introduction 



Two major bottlenecks hamper the straightforward (but toilsome) implementation of QCD 
perturbation theory in physical computations. The first that one encounters in producing 
new QCD calculations for a certain process is the evaluation of the relevant matrix elements. 
The second regards the use of those matrix elements in the actual computation of physical 
quantities. Recent years have witnessed much progress in perturbative QCD calculations 
and both these major bottlenecks have been greatly reduced. On the one hand, new 
techniques have been developed to evaluate QCD amplitudes at one loop U and on the 
other, completely general, process-independent methods 0-0] have been set up to compute 
physical cross sections at next-to- leading order (NLO). 

In particular, reasonably compact expressions for all five-parton one-loop amplitudes 



have been obtained PHTI1 and NLO calculations for physical processes involving many par- 



tons have become feasible. The three-jet cross section in hadron collisions in the simplified 



case of pure-gluon subprocesses is now available |TT|,[T2|, as is the four-jet cross section in 
electron-positron annihilation in the simplified cases of leading colour |13| and four-quark 
final states Hi ; the full QCD results are expected to appear soon. 



These two bottlenecks are, however, strongly related, and overcoming them indepen- 
dently is not sufficient for a successful implementation of perturbative QCD. The relation- 
ship is due to the regularization of unphysical divergences and, ultimately, to unitarity. 

Theoretical evaluation of jet observables, i.e. infrared and collinear safe quantities, 
should lead to unambiguous results for physical cross sections. As a statement of prin- 
ciple this is almost trivial, but its practical realization is far from trivial. This is because 
expressions at intermediate steps of the calculation (loop matrix elements and integrals of 
tree- level matrix elements) contain ultraviolet, soft and collinear divergences. Consequently, 
one has to introduce some regularization procedure. Eventually, the regularized singular- 
ities cancel in physical quantities and the finite remainder should be independent of the 
regularization procedure. This independence is only achieved, however, if the regularization 
prescription is unitary. The evaluation of the matrix elements and their implementation 
in the calculation of physical observables must therefore be carried out without violating 
unitarity. 

A regularization scheme (RS) that fulfils the requirement of unitarity in QCD is known: 
conventional dimensional regularization. However, the intermediate ingredients that are 
necessary for a complete evaluation of cross sections are not always available within this 
RS. In fact, the new techniques for calculating one-loop QCD amplitudes use a different 
version of dimensional regularization, dimensional reduction |15[ . Moreover, even when the 
squared matrix elements in conventional dimensional regularization (or any other unitary 
RS) are available, calculations in other schemes may be simpler. Thus it is worth while 
investigating whether the RS dependence can be controlled in a simple way. 

Within the context of dimensional regularization, the RS dependence of one-loop am- 
plitudes has been studied in detail. In Ref. fll6fl , from the general structure of the squared 
matrix elements for all 2 — ► 2 QCD subprocesses, effective transition rules were derived 
to relate the NLO loop corrections in conventional dimensional regularization to those in 
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several different dimensional regularization schemes, including dimensional reduction. For 
practical purposes, this information on the matrix elements is sufficient for most of the 
QCD computations at NLO. 

Using a mass regularization scheme, many subtleties related to a consistent unitary 
implementation of QCD amplitudes in cross-section calculations have been pointed out in 
Ref. |17| for the specific case of three-jet production in e + e~ annihilation. The treatment 



of more complicated processes along the lines of Ref. |17] is far from trivial. 



The general algorithm for computing jet cross sections presented in Ref. |J provides a 
very simple and transparent way to study the RS independence and unitarity of QCD cross 
sections. In this paper we explore and illustrate some of the issues. Our main results are: 

• we show how the several RS-dependent ingredients can explicitly be isolated in simple 
contributions to the NLO cross sections; these contributions are universal (they do 
not depend on either the process or the jet quantity) and provide an explicit control 
on the RS independence of the calculation; 

• in the general context of dimensional regularization, we derive the explicit RS de- 
pendence of the one-loop QCD amplitudes in a simple way without doing any loop 



calculations; in particular, we confirm the transition rules obtained in Ref. ||16|| : 

• more generally, we can provide an explicit and simple recipe to guarantee unitarity 
of cross-section calculations for any regularization prescription that is consistently 
defined at the level of one-loop matrix elements. 

As for the last point, we should point out that the RS issue considered in this paper 
does not regard the ultraviolet behaviour. Ultraviolet divergences have to be properly reg- 
ularized and then renormalized in off-shell Green functions. This leads to the introduction 
of the running coupling 0:3 (/i 2 ), which we always assumed to be defined in a fixed renormal- 
ization scheme, say, the MS scheme (the renormalization-scheme dependence can always be 
controlled by an overall perturbative shift in as). The on-shell limit of the Green functions 
thus defines the singular (because of soft and collinear divergences) loop amplitudes we are 
concerned with. These singularities are unphysical in the sense that they disappear in any 
jet observable, i.e. in any physical quantity that is well-defined in QCD perturbation theory. 
Thus there is no need to consider only (soft and collinear) regularization prescriptions that 
are completely justified on a field theoretical basis (for instance, to any perturbative order 
or for both tree- level and loop amplitudes). One can introduce any regularization prescrip- 
tion that is well-defined at the sole level of loop amplitudes: by explicit construction, we 
shall show how one can then enforce unitarity in the cross-section calculation at NLO. 

The outline of the paper is as follows. In Sect. 2 we recall the general and precise 
definitions of the NLO QCD cross sections we aim to calculate, and discuss their relationship 
with the unitarity condition. In Sect. 3 we recall the important features of a general 
formalism — the dipole formalism M - to calculate those cross sections. We identify 
the terms that may contain RS dependence, and write down the explicit conditions that 
ensure unitarity of physical cross sections for any RS that is consistently defined at the 
level of one-loop matrix elements. In Sect. 4 we show that, within the general framework 
of dimensional regularization, the explicit RS dependence of NLO QCD calculation can be 
derived from that of corresponding the Altarelli-Parisi splitting functions. We do that for 
the main schemes in current usage. Sect. 5 contains our conclusions. 
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2 QCD cross sections at next-to-leading order 



According to the QCD factorization theorem, the hadron-level cross section for a jet ob- 
servable has the following expression 

(Xhad = / x a x d , (1) 

where a is the parton-level cross section and / and d are the non-perturbative parton 
densities and parton fragmentation functions of the incoming and observed (in the final 
state) hadrons, respectively. The notation in Eq. ([]]) is symbolic (see, for instance, Sect. 6 
of Ref. for a detailed notation). The hadronic cross section <7h a d depends on the definition 
of the jet quantity and on the momenta of initial-state and final-state observed particles. 
The parton-level cross section depends on the jet definition, on the parton momenta and on 
the parton flavours. The crosses in Eq. (P stand for the convolution over the momentum 
fractions and for the sum over the flavours. 

We recall that o"h a d is a physical cross section while /, d and a are not separately 
physical quantities. They depend on the factorization scheme. Having defined this scheme, 
i.e. the process-independent parton distributions / and d, the partonic cross section o is 
computable with no ambiguities in perturbation theory to any order in as- Its perturbative 
expansion up to NLO is the following: 

^LO | „NLO (n\ 

a = a + o . [2) 

The leading order (LO) contribution a LO is obtained by integrating the fully exclusive 
cross section da B in the Born approximation over the phase space for the corresponding 
jet quantity. Suppose that this LO calculation involves m partons in the final state. Thus, 
we write 

M> - 1 da B . (3) 



Note that the LO cross section in Eq. @ is finite by definition. 

Using analogous notation, the NLO cross section a NLO is a sum of three integrals: 



a NLO = da H + da v + da c . (4) 

Jm+l Jra Jra 

Here da c is a counterterm that defines the factorization scheme, da R (the 'real' cross 
section) is the exclusive cross section with m + 1 partons in the final state, and da v (the 
'virtual' cross section) is the one-loop correction to the process with m final-state partons. 

Strictly speaking, the expression on the right-hand side of Eq. (f|) has only a formal 
meaning because its contributions are separately divergent. The virtual cross section da v 
is proportional to the one- loop matrix element that, although renormalized, still contains 
soft and collinear singularities coming from the loop integration. In order to remove these 
infinities, the loop integral has to be regularized and, correspondingly, da v is replaced by 
its (RS- dependent) regularized version da^ s . The real cross section da R is finite but its 
integration over the m + 1-parton phase space produces soft and collinear divergences that 
cancel those in da v , thus leading to a finite NLO cross sectionQcr ArL °. This finite remainder 

^The integral of da R also produces additional collinear divergences that cancel those in da c . Although 
these are in fact regularized by the same scheme as the others, the resulting RS dependence can be absorbed 
into the factorization scheme dependence of da c , which does not concern us here. 
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is unambiguously denned provided that, before its integration, the real cross section da R is 
replaced by a consistently regularized version da R & . A more correct way of writing Eq. (^) 
is thus the following 

a NLO = f da* & + / da% & + f da c , (5) 

where some formal limit R.S. — > is understood. 

Consistency of the RS means unitarity. More precisely, unitarity of the regularized 
theory implies that the scattering amplitudes T^ b & that are used to compute its matrix 
elements have to fulfil the following unitarity condition 

2ImT- = £|Tr| 2 , (6) 

b 

up to the relevant order in perturbation theory. In particular, at NLO the discontinuity 
of the one-loop matrix element (which is used in da^ & ) on the left-hand side of Eq. (f|) 
provides a constraint on the squares of the tree-level matrix elements (which are used in 
da R & ) on the right-hand side. 

The unitarity condition in Eq. @, which ultimately is at the basis of the cancellation 
theorems ||18||, also shows that the main difficulty in controlling the cancellation of diver- 



gences in Eq. (f|) has a kinematic origin. If we had to compute the total cross section 
for a given process, the real and virtual contributions could be combined at the integrand 
level before doing the loop integral. In the soft and collinear regions, this integrand would 
contain exactly the difference between the left- and right-hand sides of Eq. (P) and, hence, 
would be integrable even without introducing any regularization prescription. On the con- 
trary, in Eq. @ the integrations of the real and virtual contributions have to be performed 
over different phase-space regions. The shape of the two phase-space regions depends in 
a non-trivial way both on the number of partons and on the actual definition of the cross 
section. In principle, one has to do a detailed calculation for any different jet observable in 
any given process. 

A general strategy [0-0] to overcome this difficulty consists in trying to expose the 
cancellation of singularities directly at the integrand level. This amounts to recasting 
Eq. (|j) in the following form: 

a NLO __ a NLO{m+l} _|_ ^NLO {m} _|_ ^NLO {m} /y\ 

In the contribution a NLO { m + 1 } ; the integration is carried out over the m + 1-parton phase 
space. The contributions a NLO { m } and a- NL °^ m \ instead, involve the integration over the 
m-parton phase space. 

The precise definition of these three terms depends on the detailed method used in going 
from Eq. (|4]) to Eq. (|7]). The common feature is that the three integrands are separately 
finite. More precisely, the integrand in a NL °{ m + 1 } [ s explicitly RS independent in the sense 
that it can be defined without introducing any regularization scheme. The integrand in 
^NLO{m} j g a j go eX pii c itly RS independent, although it depends on the factorization scheme. 
The integrand in a NLO ^ contains the sum of two terms: they are separately divergent if 
the regularization is removed. Both the divergences and the RS dependence (should) cancel 
in the sum. In the rest of this paper, we concentrate on the finiteness and RS dependence 
of the integrand in cr iVLO { m }. 
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3 Dipole formalism and the issue of unitarity 



The dipole factorization formulae and the general algorithm for computing QCD cross sec- 
tions presented in Refs. fHH are particularly convenient for studying the RS dependence of 
the integrand in a NLO . They indeed provide completely general and explicit expressions 
for the cross-section contributions on the right-hand side of Eq. (|7|). In particular, all the 
kinematic complications related to the actual definition of the jet quantity are confined in 
overall factors. 

The key point of the dipole formalism is the universal definition of a 'fake' cross section 
da^s. t n& t depends on the momenta of the m + 1 real partons involved in the evaluation of 
the NLO cross section in Eq. (|j). The fake cross section has the following general form 

d°L = E ® ^pole . (8) 

dipoles 

In Eq. (H) the only dependence on the physical quantity we are interested in is contained 
in the factor da^ & . This is exactly the (regularized version of) Born-level cross section that 
enters in the calculation of the LO cross section a LO in Eq. (|3|). The only other ingredients 
needed to construct da^ & are the dipole factors dV^ olc . They depend on the RS but are 
otherwise universal: the dipole factors are completely process- and observable-independent 
and can be given once and for all || starting from the regularized S-matrix of the theory. 
The symbol <8> denotes properly defined correlations between colours and helicities of the 
partons in da^ a and in the dipole factors (see Sect. |4.1|) . 



The first main property of Eq. (||) is that, in the soft and collinear regions, da^ & has the 
same pointwise singular behaviour as the (regularized) real cross section da* & in Eq. (|J): 

da* — ► da^ . (9) 

/ soft and/or \ 
\ collinear / 



Equation (^) follows from the factorizing properties |19| of soft and collinear radiation in 
gauge theories and the dipole factorization theorem introduced in Ref. ||. The dipole 
factors are precisely defined || by the emission probability of soft and collinear partons. 

There are several dipole terms on the right-hand side of Eq. (§). Each of them corre- 
sponds to a different kinematic configuration of the m + 1 real partons. Each configuration 
can be thought of as obtained by an effective two-step process: an m-parton configuration 
is first produced and then one of these partons decays into two partons. The Born-level 
cross section da^ & depends on the m-parton configuration, and the dipole factors describe 
the one-to-two parton decays. This two-step pseudo-process can be defined without intro- 
ducing any approximation on the m + 1-parton kinematics, thus leading to the second main 
property of Eq. (H): exact factorization of the phase space ||. 

Exact factorization means that we can carry out a factorizable mapping from the m + 1- 
parton phase space to an m-parton subspace, identified by the partonic variables in da^ & , 
times a single-parton phase space, identified by the dipole partonic variables in dV^p 0lc . 
This mapping makes dV^p Qle fully integrable analytically. We can write: 

/ ^ = E / *4* ® / ^di P s oie = / ® , (io) 

Jm+l dipoles Jm Jl Jm 
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where the universal factor J Ra is defined by 



V ~ ' « ole . (11) 



/ j 1^ dipole 
dipoles 



These two main properties of da^ & allow us to obtain the decomposition in Eq. (|7|) by 
a straightforward implementation of the subtraction method |20|| . The fake cross section in 
Eq. (|) can be subtracted from da^ & and da° , and then added back to the right-hand side 
of Eq. (|). This subtraction defines the integrands of a NLO{m+i} and a NL °^ m} in Eq. (0), 
which, owing to Eq. (|9|), are explicitly RS independent, since the RS can be removed already 
in the integrand. The remaining term defines cr NL °{ m } . Indeed, combining the virtual cross 
section da^ & with Eq. (|T0D, we can write this contribution as follows: 



a NLO {m} 



da KS + da^ & 



RS=0 



(12) 



The second term in the square bracket of Eq. (|I2|) contains all the regularized singularities 
that are necessary to cancel the (equal and with opposite sign) regularized singularities 
in the virtual correction da^ a . After adding these two terms, one can thus remove the 
regularization (as implied by the notation R.S. = 0). 

Actually, as explicitly denoted on the right-hand side of Eq. (|i~2|), the regularization can 
be removed before carrying out the integration over the m-parton phase space. The reason 
for this is the following. Firstly, as recalled below Eq. (£|), soft and collinear divergences 
in da^ & arise from the loop integral in the one- loop matrix element, independently of the 
definition of the virtual cross section. Secondly, in the second term on the right-hand side of 
Eq. (|1^) the dependence on the cross-section kinematics is fully contained in dcr^ & , which, as 
pointed out below Eq. (§), is integrable by definition. Following this observation, Eq. (|l|) 
can be rewritten so as to eliminate all the kinematic dependence from our discussion. 

The virtual cross section da^ & can be written as follows: 

da** = E ^ (m) (M) l^ Ra ({p})l?-ioop , (13) 

{m} 

where J2{ m } stands for the sum over all the configurations with m partons, \-M K& ({p})\i-\ oop 
is the one-loop matrix element squared and {p} denotes its dependence on the parton mo- 
menta (these include the m final-state partons as well as possible partons in the initial 
state). The factor d^ m \{p}) in Eq. ( |13D contains all the other contributions to the dif- 



ferential cross section: spin and colour average factors, the m-parton phase space and, in 
particular, the explicit definition of the jet observable in terms of the parton momenta {p}. 

The Born-level differential cross section da^ & , being itself defined on the m-parton phase 
space, has exactly the same form as Eq. (|13D , apart from replacing the one-loop matrix ele- 



ment |-A4 RS |i_ loop with the corresponding tree-level matrix element |.M Ra | 2 . Thus, inserting 
Eq. (|13D and the analogous expression for da^ & into Eq. ([12]) , we obtain 



a NLO {m} 



[ E ^ (m) (M) H{p}) , (14) 

{m} 
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where (omitting the notation [. . .] Ra =o> from now on) we have introduced the following 
quantity 

H{p}) = l-M Ra (M)l?-ioo P + |.m rs (M)| 2 ® / Ra (M) • (is) 

The core of Eq. (0), namely the function J-"({p}), embodies all the relevant RS information. 
Although obtained by adding the two separately divergent and RS-dependent terms on the 
right-hand side of Eq. (|T5|), this function is finite and RS independent in any unitary RS. 

Equation (|15| ) is our master equation for the study of unitarity and RS independence 
of physical cross sections. It can be considered as the analogue of the unitarity condition 
(^). From a field-theory viewpoint, Eq. ( |ToD is not as basic as Eq. (|5|). Nonetheless, it 
provides us with an explicit implementation of Eq. (^j) directly at the level of cross-section 
calculations and thus, it can be quite useful in practical terms. In order to discuss this 
point, let us consider the main features of Eq. (|i"5|). 

The master function ^{{p}) is universal. It controls the soft and collinear singularities of 
physical cross sections but depends on process- and observable-independent contributions, 
the one-loop and tree-level matrix elements and the integral of the dipole factors. This 
makes explicit the statement in the first item on the list in Sect. |I[ 

Within the context of the dipole formalism and the general algorithm for computing 
QCD cross sections of Ref. 0, Eq. (|15D allows a straightforward implementation of different 
RSs in actual calculations of physical quantities. 

In general, one can use Eq. (|15D as a simple recipe to enforce and guarantee unitarity 
and RS independence in the evaluation of physical cross sections. As recalled in Sect. [I], 
progress in the computation of loop amplitudes cannot be disjoint from similar progress in 
the use of these amplitudes for cross section calculations. In principle, one should regularize 
unphysical soft and collinear divergences by introducing a RS that is manifestly unitary, 
evaluate accordingly the real and virtual cross-section contributions in Eq. (|]), and perform 
all the steps that are necessary to end up with finite physical cross sections. Equation (pf) 
can be used as a convenient short cut for this procedure. In order to actually evaluate 
\M KS ({p})\\_ Xoop1 it is sufficient to introduce a regularization prescription that, at the 
level of one-loop amplitudes, is defined in a consistent manner (it should not spoil general 
properties such as gauge invariance and renormalizability) . Then, one can compute the 
second term on the right-hand side of Eq. (|T5|) accordingly: 

1. the partons in the tree-level matrix element Ai n& ({p}) have to be treated like the 
external partons in loop amplitudes; 

2. the parent parton and its (soft and collinear) decay partons in the dipole factors 
^Khpoie h ave to be treated like the partons inside loop integrals; 

3. the dipole phase space involved in the integral (|TT| ) of the dipole factors has to be 
treated like the phase space in the loops. 

These rules are sufficient to calculate the master function T({p}) that contains all the 
relevant dynamical information on soft and collinear divergences. Having ^{{p}) to hand, 
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one can easily perform RS transformations independently of the use of the dipole formalism 
in actual cross-section computations. 

The reason why the unitarization recipe just outlined works is simple. The integral 
I K& ({p}) of the dipole factors is proportional (but with opposite sign) to the discontinuity 
of the one-loop amplitude in the soft and collinear regions. Therefore, any unitarity defect 
in the regularization of the loop integral is automatically corrected by a corresponding 
contribution in the dipole factors according to Eq. @. Note that, in this respect, the concise 
notation in Eq. (|T5| ) may appear confusing. In fact, the tree-level amplitudes Ai n& ({p}) 
are matrix elements of the regularized theory but do not necessarily provide a complete set 
of them. The matrix elements Ai K& ({p}) are evaluated over the customary partonic states, 
while the saturation of the unitarity condition (||) requires the sum over all possible states 
that contribute to the discontinuity of the loop amplitude. The regularization prescription 
of soft and collinear singularities in the loop can introduce unphysical statesjj]: in Eq. (|15|) 
the contribution of these is cancelled by analogous terms in the integral of the dipole factors. 

Note also that our unitarization recipe does not require the actual calculation of loops. 
It is sufficient to integrate the dipole factors, which are tree-level objects in every respect. 
This is a non-trivial computational simplification. 

The use of infrared regularization prescriptions that are not manifestly unitary may 
appear an oddity. However, there are in fact such schemes in current practice. In the case 
of ultraviolet divergences, for instance, an unequal treatment of particles inside the loop 
and external particles can easily be reconciled with unitarity. A similar unequal treatment 
in the case of soft and collinear divergences is not so harmless. In the calculation of physical 
cross sections, the infrared singularities of loop amplitudes are cancelled by corresponding 
singularities arising from the integration of tree-level matrix elements: in the latter, it is 
not so trivial to make a distinction between 'internal' and 'external' particles. This point 
can be clarified by the explicit examples considered in Sect. [|in the context of dimensional- 
regularization prescriptions. 



4 Dimensional regularization 

4.1 Conventional dimensional regularization 

The RS known as conventional dimensional regularization simultaneously regularizes ultra- 



violet p2| , |23[| and soft and collinear divergences [p4fl . It amounts to analytically continuing 
parton momenta to d = 4 — 2e space-time dimensions and to considering d — 2 helicity 
states for gluons and 2 helicity states for massless quarks. No distinction is made between 
real and virtual partons. This RS is manifestly Lorentz and gauge invariant and consistent 
with unitarity. 



J The role of these unphysical states in fulfilling the unitarity condition is somewhat analogous to that 
of the Faddeev-Popov (ultraviolet) ghosts. In some specific cases they act as a negative number of scalar 
fields, so we can call them 'infrared ghosts', although in general this analogy is only heuristic. Because of 
the way in which we treat them, they are more like Feynman's 'dopey particle' than the formal ghosts 
of Faddeev and Popov. 
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In Ref. 0, the calculations were carried out using conventional dimensional regulariza- 
tion. Here we summarize the explicit results. For this purpose, we recall some notation. 

It is useful to introduce a basis {\ci, . . . , c n ) <8> |si, • • • , s n )} in colour + helicity space in 
such a way that the tree-level matrix element with n partons can be written as follows: 

M cl >-> m -> Sn {pi,...,p n )= ((c 1 ,...,c n \^(s 1 ,...,s n n |l,...,n) , (16) 

where {ci, . . . , c n }, {si, . . . , s n } and {pi, . . . ,p n } are respectively colour indices, spin indices 
and momenta of the partons. Thus |1, . . . , n) is a vector in colour + helicity space^. 

According to this notation, the (RS-dependent) tree-level matrix element squared summed 
over colours and helicities is: 

l-M Ra (M)| 2 = RS (l,...,n|l,...,n) RS , (17) 



while the colour and spin correlations denoted by ® in Eq. (|T5|) (and everywhere throughout 
the paper) are given as follows: 

\M™({p})\ 2 ® I Ra (M) = M (l, . . . , n| / RS (M) |1, . . . , n) RS , (18) 

where the integral J Ra of the dipole factors is a matrix in colour + helicity space and acts 
as an insertion operator on the right-hand side of Eq. (fL8|). 

In conventional dimensional regularization, although the dipole factors are helicity- 
dependent, spin correlations vanish after integration over the dipole phase space. Thus, 
the insertion operator^] I({p}, e) is diagonal in the helicity space and depends only on the 
colour charges Tj (see Sect. 3.2 of || for their detailed definition) and momenta pi of the 
partons in the tree-level matrix element Ai({p}) (I — 1, ... ,n). Its explicit expression is 

'«* e > = -Sf(lbj? *f v 'M £ t ' Tj {&)' ' (19) 

where /z is the dimensional-regularization scale and the singular (for e — * 0) function V/(e) 
depends only on the parton flavour and has the following e-expansion: 

V I {e)=T)(^ 2 -^+ ll l - + ll + K I + 0{e) . (20) 

For present purposes there is no need to recall the detailed calculations leading to Eqs. ( |19"l , p0"l) . 
It is sufficient to note that the constants 7/ and the Kj in Eq. ( ^0[) are related to the d- 
dimensional integral of the Altarelli-Parisi splitting functions consistently evaluated in 



conventional dimensional regularization (cf. Sect. 4.3). As a matter of fact, we have 



-y^j e + G( e 2 ) , (21) 



§In the case of initial-state partons, the definition of the state vector |1, . . . , n) in Eq. ( |16|) differs by a 
normalization factor (proportional to the number of colours) with respect to the definition used in Ref. || 
(cf. Eq. (3.11) in §). 

^Since we consider conventional dimensional regularization as the default case, we drop the label RS in 
this scheme. 
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where (P &(z;e)) denotes the azimuthally averaged splitting function. 

For the sake of completeness, the actual values of the constants entering into Eq. 



arc: 



T 2 9 = C A , 

l 9 = ^rC A - 2 ~T R N f , (22) 
K„ = K„ = I - - — I C, . K, = I — - — I C \ - — TrN 





, - n,, - ( - - — r /, . n., - i — - — r * - -g- , 

where C F = (iV c 2 - l)/2N e , C A = N c , T R = 1/2, N c is the number of colours and N f is the 
number of massless flavours. 



4.2 Regularization prescriptions within 
dimensional regularization 

Dimensional regularization was invented to regularize ultraviolet divergences in loop inte- 
grals of gauge theories in a gauge-invariant manner. In this respect the essential ingredient 
is the continuation of loop momenta into d ^ 4 dimensions. Having done this, one is left 
with some freedom regarding the dimensionality of the momenta of the external particles 
as well as the number of polarizations of both external and internal particles. 



The original choice of 't Hooft and Veltman p2[ was to continue the particle momenta 



and the helicities of vector particles inside loops into d ^ 4 dimensions, while keeping the 
momenta and helicities of external particles, as well as fermion helicities inside loops, in 
four dimensions. 

Another version of dimensional regularization, namely dimensional reduction, was intro- 
duced in Ref. [PJ in order to explicitly preserve supersymmetric Ward identities in gauge 
theories. This dimensional regularization prescription consists of continuing the virtual 
momenta of loop integrals into d dimensions while keeping all (internal and external) po- 
larizations in four dimensions. Although the operational definition of dimensional reduction 
is not complete in some higher- loop computations [p5|] , it has been explicitly checked [p6 



that gauge and supersymmetry invariance is maintained up to two-loop order. In practical 
one-loop calculations this prescription works like conventional dimensional regularization 
apart from a subtle (and important) point: one has to distinguish between four-dimensional 
metric tensors coming from the Lorentz algebra of the spin indices and (i-dimensional met- 
ric tensors arising from momentum integrals with more than one loop momentum in the 
numerator. As for the momenta of the external partons one can keep them in d dimensions 
(as originally proposed and used, for instance, in Ref. |27|| ) or in four dimensions |l(|. The 
latter option, in particular, is systematically used in the string-theory inspired techniques 
|28| for computing QCD one-loop amplitudes [|l|,|8| — pTOfl : it leads to extreme simplifications 



of unnecessarily cumbersome expressions arising, for instance, in conventional dimensional 
regularization. 

All these various dimensional regularization prescriptions work as ultraviolet regulators. 
As for gauge invariance, to our understanding the key point is that, in the procedure of 
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analytic continuation, there are only two possibilites for choosing the number of external 
polarizations. One can set them either to their four- dimensional value or to the number of 
polarizations inside the loop. This prevents the propagation of spurious (additional) degrees 
of freedom from the loop to external states. When we regularize ultraviolet singularities, 
these different prescriptions simply lead to different renormalization factors in off-shell 
Green functions and, possibly, to a perturbative shift in the definition of the renormalized 
running coupling. 

As soon as the same prescriptions are used to regularize also soft and collinear diver- 
gences, unitarity has to be carefully considered in the evaluation of the on-shell matrix 
elements, as first pointed out in Ref. JTB| . 



In summary, we are going to discuss dimensional regularization prescriptions that at the 
level of one-loop amplitudes are defined as follows. The parton momenta in the loop are 
d-dimensional, while the external momenta are either <i-dimensional (conventional dimen- 
sional regularization and dimensional reduction as in P7| ) or four-dimensional ('t Hooft 
and Veltman and dimensional reduction as in fl]i~6|). Correspondingly, the number n s (g) 



of gluon polarizations in external states is either d — 2 or 2. The number h g of gluon helicity 
states in the loop is analytically continued to d dimensions, as in conventional dimensional 
regularization and in the 't Hooft- Veltman prescription: 

^cdr = h w = d _ 2 = 2 _ 2e ; (23) 

or kept fixed in four dimensions, as in dimensional reduction: 

/i° R = 2 . (24) 

The number h q of massless-quark polarization states in the loop is 2 (as in all the reg- 
ularization prescriptions discussed above) or arbitrarily continued to d dimensions by 
2h q = Trl = 4 + 0(e), Tr 1 being the dimensionality of the spinor space (the defini- 
tion Tr 1 = 4 — 4e was used in Ref. [p9j] ). This freedom in defining h q does not simplify 
any practical calculations and is correctly referred to as harmless in any textbook that 
introduces dimensional continuation as an ultraviolet regularization. In order to emphasize 
in the simplest way the differences between ultraviolet and infrared regulators, however, 
we also consider a toy scheme that is identical to conventional dimensional regularization, 
but witr|J] h q = 2 — 2e. We will see that this leads to effects that, although trivial to keep 
track of, are not harmless. 

We summarize these definitions in Table [I]. 

The implementation of these prescriptions in the master formula QT5D is straightforward. 
Following the unitarization recipe in Sect. |3|, the tree-level matrix element M n& ({p}) has 
to be evaluated considering the partons {p} like the external ones in the loop amplitude. 
The dipole phase space is always (^-dimensional. The dipole factors are obtained from 
the emission probabilities of soft and collinear partons. However, soft emission, being 
gauge invariant and independent of the spin of the radiating parton, is insensitive to the 
treatment of the spin polarizations. Eventually, in order to relate different dimensional 

"it might seem that this is the most natural scheme, since it corresponds to a theory in which every 
quantity is analytically continued to the same number of dimensions, d. 
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Conventional 


't Hooft- 


Urigmal 


Modern 


Toy 




dimensional 


Volt man 


dimensional 


dimensional 


scheme 




regularization 




reduction 


reduction 






V 1 7 I \ 


XI V 






toy 


Number of internal dimensions 


d 


d 


d 


d 


d 


Number of external dimensions 


d 


4 


d 


4 


d 


Number of internal gluons, h g 
Number of external gluons, n s (g) 


d-2 


d-2 


2 


2 


d-2 


d-2 


2 


2 


2 


d-2 


Number of internal quarks, h q 
Number of external quarks, n s (q) 


2 


2 


2 


2 


d-2 


2 


2 


2 


2 


d-2 



Table 1: Definitions of various regularization prescriptions of one-loop amplitudes referred 
to in the text. 

regularization prescriptions, we simply have to compute the corresponding Altarelli-Parisi 
splitting functions for collinear emission. 



4.3 Altarelli— Parisi splitting functions in various dimensional 
regularization schemes 



Let us consider the (time-like) splitting of a massless parton with flavour a into two massless 
partons with flavours b and c and momenta pb and p c . The collinear limit k± — > is precisely 
defined by introducing the following parametrization of the parton momenta 



Pb 



zpv + k^ 



z 2p ■ n 



(1 - z)p" - fc£ 



k 2 



1 — zip ■ n 



(25) 



where the light-like (p 2 = 0) vector p M denotes the collinear direction and n M is an auxiliary 
light-like vector that is necessary to specify the transverse component k± {k\ < 0, k±p = 
k±n = 0) or, equivalently, how the collinear direction is approached. 

The probability of the splitting process a — >• b + c (summed over colours and spins 
of b and c) in the collinear limit is proportional to the Altarelli-Parisi splitting function 
P^ b & (z, k±). Since the collinear partons b and c have to be treated like the partons inside 
loop amplitudes, this splitting function is RS dependent. Using dimensional regularization 
the momenta in Eq. (p5|) are always rf-dimensional, but the number of polarization states of 
the partons b and c still depends on the detailed regularization prescription. Moreover, the 
function P^ b & (z, k±) depends not only on the longitudinal-momentum fraction z involved 
in the splitting process fl25D but also on the transverse momentum k±_ and on the helicity 
of the parent parton a (this parton is treated like the external partons in loop amplitudes). 
More precisely, P^f is a matrix acting on the spin indices of the parton a. The calculation 



of the polarized splitting functions is straightforward p0[. We find: 



8\Pg(z,k ± )W = (s\P*f(z,k ± )\s') = 5 S8 ,C 



{sK 9 s (z,k ± )\s') 



2z 



2 9 



(1-z) 



(*\P t 



RS. 

99 



[z,k ± )\s') = 5 SS >C F 



2(1 -z) 



(26) 



(27) 
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(li\P™(z,ki_)\u) = (ji\P* & (z,k±)\u) =Try 



-g^ + 4z(l - z 



( f i\P™(z,k ± )\v) = 2C A 



-g 



l-z 



z 



rv rv 



(28) 
(29) 



where the spin indices of the parent parton a have been denoted by s, s' if a is a fermion 
and by the Lorentz indices ^, v if a is a gluon. 

The scheme dependence related to the definition of h g only affects the Altarelli-Parisi 
function in Eq. ( |28|) for the splitting process g — > g + g. For all the other prescriptions 
discussed in Sect. \i.2\ , the scheme dependence is entirely parametrized by the number of 
gluon helicity states in the loop. This dependence consistently vanishes for e — ► 0, i.e. if the 



regularization is removed. In particular, because of Eq. (23), we can immediately conclude 
that for the 't Hooft-Veltman scheme the insertion operator I K& ({p},e) entering into the 
master formula ([15]) exactly coincides with that in Eq. ([19]) for the case of conventional 
dimensional regularization. 

Owing to the helicity conservation in the quark-gluon vector coupling, the quark split- 
ting functions in Eqs. ( p6[p7[ ) are diagonal in the spin indices. The gluon splitting functions 
in Eqs. (|28| , |29[ ) turn out to be diagonal after integration over the dipole phase space. As 
a matter of fact, this integration involves the azimuthal average of the d — 2 transverse 
components, which gives**: 



rv rv | 



d-2 



pn 



(30) 



Inserting this expression into the tree- level matrix element as in Eq. (|18|), the longitudinal 
terms proportional to p^ and p v give vanishing contributions because of the Ward identity 
(gauge invariance) p^M^f" = and only the spin-diagonal term — g^ u survives. 

In conclusion, in any dimensional-regularization scheme the insertion operator I KS ({p}, e) 
is diagonal in the helicity space and has the same expression as the insertion operator for 
conventional dimensional regularization in Eq. (|P9[): 



I* a (M,e) 



Q's 



27rr(l 



E ^ v? 




(31) 



The RS dependence is embodied in the e-expansion of the flavour functions Vf & (e), whose 
coefficients are obtained by the <i-dimensional integration of the azimuthally averaged split- 
ting functions (P^ b & (z; e)): 



_?£_ + I/£S(l_ JZ ) 
1 - z 2 9 K } 



(32) 



**In the case of conventional dimensional regularization the azimuthal average coincides with the average 
over the polarizations of the parent gluon. In other dimensional regularization prescriptions, however, the 
two averages are different, due to the mismatch between the number, d — 2, of transverse components and 
the number of helicity states of the external gluons. 
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19 



2(1^) + 1_ 

z 2 9 



(P^(z;e)) = (P^(z;e))=T i 



fl'l 



1 - 



d- 2 



(P-(z;e)) = 2C, 



l-z 



+ 



d - 2 



z(l-z) 



(33) 
(34) 
(35) 



Equations (32]-|3~5) differ by terms of order e with respect to their well-known versions 
in conventional dimensional regularization (h g = d — 2, h q = 2). According to the 
notation in Sect. |4.1| , we drop the label R.S. in conventional dimensional regularization 
and parametrize the deviation from this result by the residual scheme- dependent function 
Vfb (z) as follows: 



(P^(z\ c)> = (P ab (z; c)> + 26 P*f (z) + 0(e 2 ) . 



(36) 



Obviously, all the functions V^ K (z) as well as V^(z) vanish. 

In the toy scheme, which differs from conventional dimensional regularization simply 
because of the dimensionality of the gamma matrices in quark loops, the only non- vanishing 
P-functions are: 



Pl7(z)=Vlf(z) 



2 Tr 



z 2 + (l- zf 



(37) 



For dimensional reduction we have: 

V™(z) = Vfl(z) = C F 1 -^, V™{z) = V™(z) = C F Z -, (3* 



K( z ) = C a zQ- ~ *), PSTiz) = PS?{z) = • (39) 

In this case, some other comments are in order. Since h® R = 2, the quark splitting func- 
tions (P° R (z;e)) and (P™(z;e)) in Eqs. ( j3^ , |33| ) are actually independent of e, while the 
gluon splitting functions (P™(z; e)) and (P™(z; e)) are not. In spite of having used only 
four-dimensional polarizations, the e-dependence enters into the gluon splitting functions 
through azimuthal correlations, whose average has to be carried out in d — 2 transverse 
dimensions. 

Note also that the difference between the splitting functions in conventional dimen- 
sional regularization and those in dimensional reduction can easily be understood on field- 
theoretical basis. As discussed in detail in Ref. ||26|| , the relationship with conventional 
dimensional regularization can be investigated by splitting the four-dimensional gauge field 
of the dimensionally- reduced Lagrangian £ DR into a d-dimensional vector plus 4— d = It 
additional components. This leads to the decomposition £ DR = £ CDR -f where £ CDR 
is the term that involves only the (i-dimensional vector field and coincides with the custom- 
ary Lagrangian in conventional dimensional regularization. The additional 2e components 
enter into and behave as scalar gluons interacting with the fermion and vector fields. 
The ultraviolet role of these 2e scalars was pointed out in Ref. |[26|| . For instance, the exact 
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relation between the customary MS coupling and the minimally-subtracted renormalized 
coupling in dimensional reduction can be derived by a simple calculation as the effect of 
the additional 2e scalars. Equation (|36|) displays the infrared role of these 2e 'ghost states'. 
Indeed, the functions Pj R (», V® R (z) and V gg R (z) in Eqs. ©0) are exactly the Altarelli- 
Parisi probabilities for the splitting processes q — > q(z) + 0(1 — z), q — > <t>(z) + q(l — z) and 
g — > <p(z) + 0(1 — z), where is a scalar gluon. 

Finally, we recall that setting Cp = Tr = Ca in a tree-level QCD calculation leads to 
recovering the results in N = 1 supersymmetric Yang-Mills theory pifl . The quark is re- 
placed by the gluino g and the corresponding four- dimensional Altarelli-Parisi probabilities 
fulfil a well-known supersymmetric Ward identity, namely P gg (z) + P gg (z) = P gg (z) + P gg (z) . 
As for the <i-dimensional Altarelli-Parisi splitting functions in Eqs. (j32}-|35|), their version in 
conventional dimensional regularization violates a similar Ward identity. Using dimensional 
reduction, instead, we have 

(P™& e)) + (P^(z; e)> = (Pjf (z; e)> + (P^z; e)> , (40) 

thus leading to equal decay probabilities for the two supersymmetric partners in any num- 
ber d of space-time dimensions. This result is expected for a supersymmetric regularization. 

We note that in the toy scheme, in which the number of quark states is analytically 
continued to be the same as the number of gluon states, h q = h g = 2 — 2e, the same Ward 
identity is also recovered: 

e )> + <3r (*; e )> = <4 oy 0*; e )> + ( p m& e )) ■ ( 41 ) 

4.4 The master function T and RS independence 
of the cross sections 

The terms of order e, which arise in the splitting functions computed with different regu- 
larization prescriptions, combine with 1/e-poles coming from collinear singularities in the 
integrals of the dipole factors and thus provide the insertion operator I K& in Eq. (|I~5"D with 
finite corrections. 

More precisely, as discussed in the previous subsection, the RS dependence of the in- 
sertion operator J Ra in dimensional regularization is entirely taken into account by the 
flavour functions Vf a (e). The coefficients of the e-expansion of the functions Vf & (e) are 
related to <i-dimensional integrals of the splitting functions, as in Eq. (EH|). Different regu- 
larization prescriptions lead to contributions of the order of e to the azimuthally-averaged 
splitting functions (P^ a (z;e)) and, in turn, to Eq. fl2"ip. These contributions produce a 
scheme-dependent shift in the constants K a on the right-hand side of Eq. ( pi]) and, hence, 
non-singular corrections in Vf s (e). Parametrizing these corrections as differences with 
respect to the flavour functions Vj(e) in conventional dimensional regularization, we have 

VT & (c)=V I (c)-7r + 0(c), (42) 

where the scheme-dependent coefficients 7f s are obtained by inserting Eq. (|36|) into Eq. (|2lD : 

7f = Ef^») • (43) 
b Jo 
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Using Eq. (f43l) one can straightforwardly evaluate 7f s in different RSs. In particular, 
we have 

7? V = , (44) 
and, using Eq. ( |37|) and Eqs. ( |38l , |39~D , we respectively find 

^toy = -toy = q ? lT = ~\ T R N i > ( 45 ) 

~DR ~DR 1 r< ~, DR f/K^ 

lq =lq = 2 F > 1 9 = 6 ^ ' 



Equation ( f4"2"|) can be inserted in Eq. ( pTTD to derive a simple expression for the RS de- 
pendence of the operator J RS in dimensional regularization. As a matter of fact, neglecting 
(9(e) corrections and using colour-charge conservation, J2j^iTj = —Tj (the insertion op- 
erator acts onto the vector |1, . . . , n) Ra , which is a colour-singlet state and, hence, fulfils 
the property J2jTj\1, . . . ,n) KS = 0), we obtain 

/ RS (M,e) = /(M, e )-gE7f a + W • (47) 



In terms of the master function in Eq. dig) , we can write: 

as 
2tt 



H{P\) = \M™{{p})\U oop \M™({p})\ 2 £7? s + \M™{{p})\ 2 ® I{{p},e) . (48) 



As usual, !{{p}, e) refers to conventional dimensional regularization and is given in Eq. ([191) . 



Let us first comment on the result in Eq. (fl8"D in the context of the toy scheme in 
which h q = 2 — 2e. The renormalization of the ultraviolet divergences in the off-shell loop 
amplitudes is performed exactly as in conventional dimensional regularization. The only 
difference regards the definition of the running coupling. Using modified minimal sub- 
traction, one introduces a renormalized coupling that is related to the customary MS 
coupling by <% y = as[l — a^TuNf/^K + C?(a§)]. This is the harmless dependence on the 
ultraviolet regularization: the only physical consequence in an overall redefinition of the 
coupling constant. As soon as we consider the renormalized on-shell one- loop amplitude 



on the right-hand side of Eq. (|48), we encounter soft and collinear divergences that are 
cancelled by the insertion operator I({p},e), computed in conventional dimensional regu- 
larization. However, according to our calculation, the remaining contribution has still to be 
corrected by the finite 7) oy -terms in order to give a function T consistent with unitarity and, 
hence, RS-independent cross sections. This additional unitarity correction does depend on 
the process through the tree-level matrix element and its flavour topology (7* oy 7^ 7g° y )- 

In the context of the 't Hooft-Veltman and dimensional-reduction prescriptions, the 
scheme dependence of the one-loop QCD amplitudes was studied in Ref. [TO. On the 
basis of the explicit evaluation of the one-loop corrections to all 2 — >• 2 QCD subprocesses, 
transition rules to relate these schemes to conventional dimensional regularization were 



derived (and argued to be universal). They were confirmed by the calculation in Ref. [[LO 
We can use these rules as a consistency check of our calculation and, in general, of our 
unitarization recipe. When inserted in Eq. (f£§), the coefficients tP r in Eq. ( |46| ) provide 
the master function T with the contributions that are necessary to exactly cancel the RS 



dependence of the one-loop amplitudes |A / t RS ({p})|i_ loop as computed in Ref. |TB 
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We can turn this argument round. Writing Eq. fl48|) as follows: 



l-M Ra (M)lLioo P = -l-M Ra (M)l 2 ®/(M^) + gl-M RS (W)r E7? S +-^(M) , (49) 

and knowing that the function T is RS independent, we can compute the scheme depen- 
dence of the one-loop amplitudes |-A / ( Ra ({p})|i_ loop from our calculation of the coefficients 



7f a . This is a way to rederive the transition rules found in Ref. [|16|]. Since the 7f a 's 



are obtained from the integral of the dipole factors, this derivation is universal (the dipole 
factors do not depend on any specific QCD amplitude) and essentially involve a tree-level 
calculation (the emitted partons in the dipole are on shell) instead of a loop one. 

Concluding this section on dimensional regularization, we should add a marginal com- 
ment. The RS dependence of the Altarelli-Parisi splitting functions may, in principle, 
also affect the factorization-scheme dependence of the NLO partonic cross sections if the 
collinear counterterm in Eq. (£|) is not defined accordingly (for instance, if, independently 
of the regularization procedure, one uses the same minimally-subtracted expression for 
dac)- In the notation of Eq. (|7|), this dependence is embodied in the contribution (j NL0 ^ m ^ 
and, in particular, within the dipole algorithm, it is explicitly taken into account by the 
colour-charge insertion operators K and H of Ref. ||. If one is interested in changing 
factorization scheme, jj NL °{ m } (e.g. the operators K and H) and the parton distributions 
/ and d in Eq. ([!]) have to be varied consistently. We think that pursuing the study of 
different ci-dimensional factorization prescriptions has no particular relevance, in practice. 



5 Conclusion 



In this paper we have considered some unitarity issues related to the regularization of 
unphysical soft and collinear divergences in perturbative QCD computations. Our analysis 
has been performed in the framework of the dipole formalism. It allows one to work out 
a general discussion while providing an explicit implementation of unitarity constraints in 
the calculation of QCD cross sections at NLO. 

In Sect. |3| we have shown that the regularization-scheme independence of physical ob- 
servables is controlled by the following master function 

T{{p}) = |Al Ra (M)|Lioop + \M™{{p})\ 2 ® I K& ({p}) , (50) 

which involves the one-loop and tree-level matrix elements and the integral I KS ({p}) of the 
dipole factors. The latter are obtained by the factorization formulae for soft and collinear 
emission. Using a unitary regularization scheme, infrared singularities as well as finite 
scheme-dependent remainders consistently cancel by combining the several contributions 
on the right-hand side of Eq. flffflp . 

This result can also be used to relate one-loop amplitudes in different regularization 
schemes without explicitly carrying out any loop calculations. Turning Eq. (|50| ) round, 

l-M RS (W)lLioop = -l-M Ra (W)l 2 ®/ Ra (M) +H{p}) , (si) 
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and exploiting the scheme independence of T{{p)), one can obtain \ Ai KS ({p}) |?_ loop in a 
different scheme by simply evaluating the difference of the corresponding dipole factors. 



More generally, our explicit construction of Eq. (p0|) allows one to compute physical cross 
sections also using regularization prescriptions that are not manifestly unitary. As long as 
the regularization procedure is consistently defined at the level of one-loop amplitudes, 
we can give a recipe (see Sect. to calculate the dipole-factor contributions that are 
necessary to guarantee unitarity. From the regularized virtual corrections one can thus 
extract the essential physical information, that is, the regularization-scheme independent 
function jF({p}). 

In order to make our general discussion more definite we have considered, in Sect. f|, 
the case of dimensional-regularization prescriptions. This analysis required the explicit 
calculation of the rf-dimensional Altarelli-Parisi splitting functions in various regulariza- 
tion schemes. Our results for these functions differ from previously published ones in the 
't Hooft-Veltman and dimensional reduction schemes PJIBfl. These differences emerge in 
the e-dependent terms as can be seen by comparing our formulae (|36^p9|) with the corre- 
sponding results in Refs. P,[16|j. The differences come about because in these references 
the average over polarizations of the parent parton was taken, instead of the azimuthal 
average. These are not the same in the case of gluon splitting functions as discussed in 
Sect. ^L^. However, this slight error has no practical consequence on the main results of 
Refs. [p[|l6f because in those references the relationship between splitting functions and 
RS dependence was not exploited (in Ref. (T^|, the splitting functions were used to relate 
the definition of different factorization schemes). On the contrary, we make full use of 
this relationship, allowing us to derive the unitarity corrections needed to relate different 
dimensional-regularization schemes (see Eq. (f49[)), without having to make any loop calcu- 
lations. Our results are in full agreement with those obtained in Ref. [IB|. We are able to 
give a probabilistic interpretation (see Eq. (|43|)) of the coefficients ^yf a 



The method can also be applied to other regularization procedures, such as, for instance, 
massive regularization of the loop integrals. Many higher-order calculations in QED have 
been carried out using this regularization scheme. Dimensional regularization is certainly 
preferred in massless QCD and its computational advantages have also been exploited in 
recent QED achievements . Our unitarization technique and the extension of the dipole 



formalism to massive partons [33| can be convenient in order to use known QED results 
in QCD applications and to combine them with new QED calculations in dimensional 
regularization. 

Systematic QCD calculations at next-to-next-to-leading order (NNLO) for jet observ- 
ables will become feasible only when efficient techniques for evaluating two-loop matrix 
elements will be set up. At that point the unitarity issue discussed in this paper will show 
up again. Our results may eventually be very useful to tackle this issue, provided the 
validity of the dipole formalism is extended to such a level of accuracy. 

In concluding, we would like to point out another feature of Eq. fl50|). The operator 
J Ra ({p}) is obtained by integrating the dipole factors and, correspondingly, |A / i RS ({p})|f_ loop 
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is the result of the integration over the loop momentum: 

FlM)=f ^ Ra (M)lLioo P + E / \M™{{p})\ 2 ®dVl* ole . (52) 

Jlo °P dipolos Jl 

It is conceivable that one may find a way of combining the two integrands such that the 
dipole factors act as a local counterterm for the loop integral. Achieving this, one could 
avoid the introduction of any soft and collinear regularization and the ensuing unitarity 
problems. Most importantly, one would be able to carry out NLO calculations of physical 
cross sections by the sole use of numerical methods, without any analytical calculation of 
one-loop amplitudes. 
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